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Numerical predictions of the Mars Science Laboratory reaction control system jets 
interacting with a Mach 10 hypersonic flow are compared to experimental nitric oxide 
planar laser-induced fluorescence data. The steady Reynolds Averaged Navier Stokes 
equations using the Baldwin-Barth one-equation turbulence model were solved using the 
OVERFLOW code. The experimental fluorescence data used for comparison consists of 
qualitative two-dimensional visualization images, qualitative reconstructed three- 
dimensional flow structures, and quantitative two-dimensional distributions of streamwise 
velocity. Through modeling of the fluorescence signal equation, computational flow images 
were produced and directly compared to the qualitative fluorescence data. Post processing of 
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the experimental and simulation data enabled the jet trajectory to be extracted for a more 
quantitative comparison. The two-dimensional velocity fields were reconstructed through 
interpolation of a series of single-component velocity profiles. Each distribution of single- 
component velocity was obtained using molecular tagging velocimetry. After validating the 
numerical model, the numerical solution was further examined to gain insight into 
hypersonic jet-in-crossflow interaction. Novel nitric oxide planar laser induced fluorescence 
experiments are proposed based on this analysis. 

Nomenclature 

AoA = vehicle angle of attack [°] 

D = capsule model diameter [mm] 

P = fluid pressure [Pa] 

Sf = fluorescence signal [arbitrary units] 

T = fluid temperature [K] 

% NQ = nitric oxide mole fraction 

I. Introduction 

r J^HE Mars Science Laboratory (MSL) vehicle was recently launched by NASA in November, 2011, and 

successfully performed entry maneuvers through the Martian atmosphere in August, 2012. Derived from the Viking 

mission design, the MSL is a lifting, capsule-type entry vehicle with a heat shield cone half angle of 70 degrees. The 

vehicle underwent an ambitious entry, descent, and landing (EDL) process, which included the execution of a series 

of bank reversal maneuvers using the vehicle’s reaction control system (RCS) jets [1]. As shown in side view of the 

MSL vehicle in Fig. 1, the RCS jet pairs are distributed and orientated on the aftbody to provide pitch, roll, and yaw 

moments to the vehicle. In the figure, the origin for the Cartesian coordinate system (x, y, z) is located on the nose of 

the forebody heat-shield. The aerodynamic maneuvers were used to improve the down-range and cross-range 

landing accuracy while the strategic firing of RCS jets counteracted aerodynamic instabilities and atmospheric 
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dispersions. A significant experimental and computational effort, motivated by the NASA Engineering and Safety 
Center (NESC), was undertaken to understand and quantify any large and unexpected aerodynamic interferences 
that could have occurred when the RCS jets were activated. The NESC study objectives were to validate the CFD 
models to show that there were no serious discrepancies in the aero/RCS interactions that had been predicted and to 
advance the CFD aero/RCS predictive capabilities. There is an expected interaction between the reaction control 
system (RCS) jets and the capsule wake in addition to interaction between individual RCS jets [1-4]. The hypersonic 
flow field consists of a strong bow shock wave ahead of the capsule forebody, a supersonic expansion fan around 
the capsule shoulder, and a large subsonic separated region in the capsule wake. Recompression shock waves 
redirect the outerflow downstream of the capsule wake. The capsule has a radial center-of-gravity offset that causes 
the vehicle trim at a predetermined angle-of-attack (AoA), providing lift for maneuvering. The RCS jets are used to 
perform bank maneuvers, changing the pointing direction of the lift vector to fly a guided entry. The under- 
expanded RCS jets interact with both attached supersonic cross-flows and separated subsonic recirculation zones 
depending on the jet location and vehicle AoA. Further details of the vehicle aerodynamics during entry are 
discussed in the literature [3,4]. 


Heatshield- 
Forebody (HF) 



Coordinate 

System 

Origin 


Heatshield-Shoulder-Junction (HSJ) 
Forward-Backshell-Surface (FBS) 
Backshell-Hip-Junction (BHJ) 

RCS Jet A1 

- Backshell-lnterface-Plate(BIP) 

Parachute-Closeout- 
Cone (PCC) 

RCS Jet B2 

Aft-Backshell-Surface (ABS) 


Fig. 1 MSL aeroshell configuration and nomenclature. 


There is a strong reliance on computational fluid dynamics (CFD) to aid in the design of capsule-type entry 

vehicles because of the high costs of ground-based tests and the impracticality of a flight test program in the Martian 

atmosphere. Although the computational resources available to simulate the high-speed aerodynamics are 

improving, there is still a strong dependence on simplified numerical models to predict the transport of mass, 
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momentum, and energy. Turbulence modeling, for example, needed to be included in the simulation of the MSL 
vehicle aerodynamics because of the high Reynolds numbers that were expected to occur during the peak heating 
loads of atmospheric entry [5]. High Reynolds numbers occurred because of MSL’s relatively large capsule 
diameter and high ballistic coefficient in comparison to previous vehicles designs [6]. Although direct numerical 
simulation (DNS) offers the most accurate means to predict this type of flow, the computational resources available 
are insufficient to resolve the entire turbulence spectrum. By partially modeling the turbulence spectrum, large eddy 
simulation (LES) offers a compromise between DNS and Reynolds Averaged Navier Stokes (RANS) simulations. 
Although LES alleviates the computational load requirements, the technique still requires high spatial and temporal 
resolution and long run times to collect statistics. By modeling the entire turbulence spectrum, RANS simulations 
can provide steady-state solutions of the mean variables with relatively short turn-around times. A recent 
computational study of the MSL RCS jet flow field was performed using detached eddy simulation (DES), which is 
a hybrid LES/RANS model [7]. In that work, LES was used to predict the unsteady separated flow in the capsule 
wake and the RANS-type model was used to predict the flow near the vehicle surface. That work provided insight 
into the unsteady nature of the RCS jet-crossflow interaction. Each configuration, however, took approximately 66 
hours to complete using each of the 240 cores of the University of Minnesota Supercomputing Institute’s SGI Altix 
XE Linux Cluster. Based on the long run-times, the technique is not yet practical for design optimization or for 
parametric studies. Therefore, the current strategy for future entry vehicle design is still to use RANS-based models 
for prediction and rely on a robust and comprehensive validation against trusted ground-based experiments [2,3,8]. 
As a result of the NESC initiative, recent detailed experimental data of flow separation and jet interaction associated 
with the MSL vehicle are now available for a rigorous CFD validation [4]. Validation is crucial, especially for 
RANS-type models, since they have difficulty in predicting turbulent flows with streamline curvature or flow 
separation [9]. 

In high-speed aerodynamic experiments, non-intrusive laser diagnostic measurements are preferred over 
conventional probe-based techniques to build a validation dataset. Laser-based measurements are preferred because 
it is often unfeasible to place a probe inside a hypersonic flow field since the mere presence of a physical object can 
artificially distort the flow and contaminate the results. Furthermore, some laser diagnostic techniques, such as 
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planar laser-induced fluorescence (PLIF), are able to extract highly resolved and spatially correlated data, which is 
preferred over traversing a single-point measurement probe. Quantitative PLIF techniques can also be used to 
extract flow variables such as velocity, temperature and species concentrations needed for direct comparison to 
CFD. A major problem that occurs with numerical model validation using quantitative laser diagnostics is that the 
experimental techniques are usually restricted to narrow operating conditions [10]. With computational flow 
imaging (CFI), both the fluid dynamics and spectroscopic process are modeled to produce simulated PLIF images, 
which can then be compared to experimental PLIF images. This strategy offers comprehensive model validation and 
simplifies the laser diagnostic technique and experimental setup. Some early work using the CFI technique to 
compare simulation predictions to experimental PLIF images for high-speed aerodynamic applications can be found 
in Refs. [11] and [12]. 

In the present work, numerical predictions of the MSL RCS jet flow field using the OVERFLOW code and the 
Baldwin-Barth one-equation turbulence model are compared to an experimental PLIF validation dataset. The 
experimental validation dataset was produced from a wind tunnel test in NASA Langley Research Center’s 31 -inch 
Mach 10 facility in 2010 [4]. In that study, nitric oxide (NO) gas was supplied to both the RCS jet plenum and to the 
sting that held the 1/2954 scale model of the MSL flight vehicle. This facilitated NO-PLIF visualization of both the 
capsule wake and the RCS jets for jet-on and jet-off conditions. A CFI model was applied to the numerical solution 
to enable comparisons to the experimental flow visualization data. Three dimensional distributions of PLIF signal 
intensity near the RCS jets were reconstructed from the original experimental two-dimensional PLIF images. 
Similarly, three-dimensional jet structures were extracted from the simulated PLIF field for comparison. Further 
post-processing and analysis of both the experimental and simulated flow visualization data have been performed to 
make comparisons of the RCS jet trajectory. To supplement flow visualization, molecular tagging velocimetry 
(MTV) was used to map the aftbody velocity field. After demonstrating the accuracy of the simulation, the 
numerical solution was further examined to gain insight into hypersonic jet-in-crossflow interaction and to further 
understand the limits of the PLIF measurement technique. Based on that analysis, future experiments using the NO 
PLIF technique are proposed. 
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II. Wind Tunnel and Model Experimental Setup 


The experimental data used for the validation was performed in NASA Langley Research Center’s 31 -Inch Mach 
10 Air Tunnel facility in 2010 [4]. By electrically heating compressed air, the facility is able to provide Mach 10 
flow through the 31 -inch square test section for approximately 90 seconds [13]. The test section has optical access 
through two large windows located on the top and bottom tunnel walls. A third smaller window is located on the 
side wall opposite to the model injection system. Each window transmits ultraviolet light, which allows the 
projection of the laser sheet through the top window into the test section. Digital images of the NO fluorescence 
were acquired through the side window. 

Figure 2 shows the side view of the 1/2954 scale titanium MSL entry capsule model. The diameter of the model 
is D = 152.4 mm. As shown, the capsule is at an AoA of -20 degrees with respect to the flow. This angle 
corresponds to a horizontal sting. A water-cooled 5 -component force and moment balance was located inside the 
capsule and aligned with the model’s center of gravity (obviating the need for axial force to perform moment 
transfers from the balance moment reference point to the capsule center of gravity). The balance cabling, 
thermocouple cabling, and N 2 /NO plumbing were protected from aerodynamic heating during the -90 second 
runtimes by a protective shroud. The sting AoA was adjusted during the run using a motorized control system. For 
most of the runs, NO was seeded into the flow only through the B2 RCS jet pair (Fig. 1). The RCS jet fluid was 
seeded with 2-5% NO, with the remaining gas being N 2 . This allowed observation of the shape, structure and 
trajectory of the RCS jets. In some runs, the RCS jets were not operated. Instead, pure (100%) NO was seeded into 
the flow through a gap between the sting and the shroud placed over the sting. As a result, the NO was entrained in 
the separated wake flow behind the capsule, allowing for the visualization of the wake’s shape and size. For the 
MTV measurement in the wake with the B2 RCS jet pair on, a 2-5% NO/N 2 mixture was seeded into both the jet 
plenum and model-sting junction. Reference [4] should be consulted for further details of the experimental setup and 
run conditions. 
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Fig. 2 Side view photograph of MSL model, sting with protective heat shroud, and tunnel side wall. 

Reproduced from Ref. [4]. 

III. PLIF Visualization and Velocity Measurement Techniques 

The NO PLIF system uses an ultraviolet laser sheet to interrogate a slice in the flow containing seeded NO. This 
UV light sheet excited fluorescence from the seeded NO molecules. The fluorescence was detected by the high- 
speed PI-MAX II intensified CCD digital camera. The laser sheet could be translated across the model, thereby 
visualizing different cross sections of the flow. In post-processing, these different spatial measurements could be 
recombined to show the three dimensional structure of the wake flow and RCS jets. Details of the measurement 
system used in this experiment, as well as the experimental procedure used, can be found in Refs. [14] and [15]. 

Streamwise velocity measurements were made by measuring the displacement of a series of lines aligned 
vertically with respect to the freestream flow in a streamwise plane. The lines were formed by passing a portion of 
the laser sheet through a 50-mm-long LaserOptik difhision- welded lens array. The lens array consisted of 25 x 1 m 
focal length cylindrical lenses that focused the sheet into 25 collimated beams. The beams were then passed into the 
flow and excited fluorescence in the NO molecules. A Cooke DiCAM-Pro dual-frame intensified camera was used 
to probe the molecules at the instant they were tagged (initial image) and 500 ns after excitation (delayed image). 
The initial and delayed images were compared to determine the streamwise (F*) component of velocity. A detailed 
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discussion of the experimental methodology can be found in Ref. [4]. A detailed discussion of the velocity and 
uncertainty analysis can be found in Ref. [16]. In this paper, modifications have been made to the methodology 
described in Ref. [16] for calculating the estimated time delay and timing uncertainty for sequential images. A 
difference of less than 1% in both measured mean velocity and uncertainty the in mean velocity is expected using 
these modifications versus using the methodology of Ref. [15]. The two-dimensional velocity fields were 
reconstructed through interpolation of the single-component velocity distributions. The flow visualization and 
velocity data were imported into Tecplot 360® for comparison with the CFD results. 

IV. Computational Fluid Dynamics Simulation 

CFD simulations of the flow field were performed with the OVERFLOW (the OVERset grid FLOW solver) 
Version 2.2 code [17,18]. The simulations solved the compressible Reynolds Averaged Navier Stokes equations. 
The Baldwin-Barth one-equation model was used for turbulence closure. The implicit time advancing scheme used 
was the Symmetric Successive Over-Relaxation (SSOR) method. The third-order HLLC (Harten, Lax, van Leer, and 
Einfeldt) numerical flux function with the Van Albada limiter was used. Three grid-sequencing levels were used to 
accelerate steady-state convergence. Solutions were considered converged when the first three significant figures of 
the component (i.e. backshell) mean aerodynamic force coefficient did not change for approximately 2000 iterations. 

The domain geometry, as shown in Fig. 3, is approximately three capsule diameters in the y - and z - directions 
and three and a half capsule diameters in the x - direction. The far- field boundaries are modeled with a characteristic 
inflow/outflow boundary condition. The gas composition at the inflow boundary is air. The inflow stagnation 
pressure, stagnation temperature, and Mach number are 9.11 MPa, 1032 K, and 9.94, respectively. The solid walls 
are modeled with a no-slip, adiabatic boundary condition with an estimated y + of one using a flat plate 
approximation. The RCS nozzle was modeled back to the plenum, where total pressure and temperature conditions 
from the test are enforced. At the plenum boundary, pure N 2 with a stagnation pressure of 689.5 kPa and stagnation 
temperature of 300 K was specified. The domain was discretized using overset structured grid blocks whose 
interfaces are interpolated via stencils defined using the Pegasus5 code. Approximately 50 million nodes were used 
in the simulation. 
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Fig. 3 Computational domain. 

The perfect gas law was used in the simulations. In the analysis, a CFD-computed map of the streamwise 
component of velocity was compared directly with experimental velocity measurements. Two additional steps were 
required to compare the PLIF flow visualization data with CFD predictions of the RCS jet plumes. First, the RCS jet 
gas was assumed to be 5% NO / 95% N 2 instead of pure N 2 . To clarify, pure N 2 was used in the simulations and NO 
was added later in post-processing by assuming a 5% molar fraction with respect to N 2 . Chemistry was assumed to 
be frozen. Second, the CFD solution was used to compute a theoretical PLIF image using the polynomial method 
described in Ref. [19]. Such computed PLIF images were compared with experimental PLIF intensity images using 
iso-contours in Tecplot. 


V. Validation Results 

Below, direct comparisons between CFD and the experimental PLIF data are made for three different types of 
measurements: separated wake flow measurements, RCS jet structure and trajectory measurements, and streamwise 
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velocity measurements. Comparison of the separated wake flow and RCS jet structure are qualitative while 
comparisons of the jet-trajectory and streamwise velocity are quantitative. 

A. Wake Flow Visualization without RCS Jets 

Figure 4a shows an NO fluorescence image of the separated wake behind an MSL model vehicle obtained in the 
previous study of the vehicle’s aerodynamics [4]. In this configuration, pure NO was seeded into the flow through a 
gap between the sting and the shroud placed over the sting. The RCS jets were not activated in these tests. The 
reversed flow in the separated wake convected the NO gas upstream of the Parachute Closeout Cone (PCC) near the 
separation point. The laser sheet interrogated the flow vertically from top to bottom and was positioned on the model 
centerline. Although both the laser position and model AoA were parameters that were varied in the experiments, all 
comparisons in this current study are evaluated at a model AoA of -20°. In the previous work the PLIF data were 
visualized using the Virtual Diagnostics Interface (ViDI) method, which allows for accurate placement of the 
processed image data relative to the MSL model. In the current study, Tecplot was used as the visualization tool. 
ViDI, however, was used to prepare the PLIF data before being imported into Tecplot. Figure 4b shows the 
separated wake at the same conditions described above, but with stream traces overlaid from the OVERFLOW CFD 
solution. Unlike ViDI, which displays a series of individual images, Tecplot displays a true three-dimensional 
reconstruction of all of the PLIF data. The image shown in Figure 4b is a slice of experimental PLIF data extracted 
from the three-dimensional reconstruction to show the original two-dimensional distribution of the PLIF signal. To 
clarify, after the two-dimensional images were loaded into an interpolated three-dimensional Tecplot data file, a 
two-dimensional slice was extracted to produce the image. Although the original two-dimensional images could be 
used to produce Fig. 3b, the process described above was implemented to verify the accuracy of the three- 
dimensional data files that were constructed in Tecplot. Both the field of view and image resolution were reduced to 
produce manageable three-dimensional Tecplot files to be loaded onto a desktop PC. As shown in this qualitative 
comparison, the location and size of primary and secondary recirculation zones agree well with streamline 
predictions from the OVERFLOW code. 


10 


American Institute of Aeronautics and Astronautics 



(a) (b) 

Fig. 4 Separated wake flow at model centerline (jet off): a) Experimental PLIF distribution, reproduced from 
Ref. [4] ; b) Experimental PLIF contour plot with overlaid streamlines predicted by OVERFLOW. 


By scanning the laser sheet in the spanwise (- y ) direction, the three-dimensionality of the wake can be observed. 
Figure 5 shows an isometric view of the MSL model with two-dimensional contour maps of the PLIF signal at three 
separate spanwise locations. A series of stream traces predicted by OVERFLOW are overlaid onto the figure. 
Similar to the flow separation that occurs from the backshell-interface-plate (BIP) near the model centerline (y = 0 
mm), NO becomes entrained into a larger recirculation zone that extends from the PCC to the heatshield-shoulder- 
junction (HSJ) near z = 0 mm. Separation does not occur at the top HSJ near the model centerline (y = 0 mm) since 
the model is pitched forward. This creates a relatively shallow angle between the freestream flow and model surface 
allowing the flow to remain attached. Discussion of the three-dimensional wake shape in relation to the effective 
surface angle is discussed in past work [4]. 



(a) (b) (c) 

Fig. 5 3D visualization of the separated wake flow through translation of laser sheet: a) y = -20 mm; b) y =-46 

mm; c) y = -63 mm. 
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Note that the experimental PLIF data do not extend far enough in the downstream (+x) or in the vertical (-z) 
directions to fully capture the size of the separated wake. Restrictions in the camera field of view and laser sheet 
dimensions prohibited these measurements. Where the experimental data are available, the OVERFLOW prediction 
of stream traces agrees well with the outer boundary of PLIF signal. Near the sting, the predicted stream traces 
indicate that the large recirculation zone forms into a trailing vortex. Although vorticity in the near field is 
dominated by the z- and y- component of vorticity, redistribution to x-component vorticity occurs in the far field. 
Although the PLIF data are not available downstream of the main separated region, NO is expected to exist further 
downstream because of entrainment into the trailing vortices. 

Figure 6 shows a three-dimensional iso-surface of the PLIF signal with the simulated stream traces. To compare 
the overall shape of the separated region, side, top, back, and isometric views are shown in Fig. 6. The vertical and 
streamwise limits of the PLIF data field-of-view are evident from the images. For example, full separation is 
predicted on the lower forward-backshell surface (FBS) and over the aft-backshell surface (ABS), but PLIF data are 
not available for comparison (denoted by the hard edges in the iso -surface data). For the data that are available, there 
is good agreement between the experimental and computed separated shape. The simulation predicts a dip in the 
separated shape that is confirmed by the experimental data. The location of the dip is indicated by the red arrow in 
the isometric view of the figure. This location corresponds to the interface that forms between the counter-rotating 
vortices in the far field. The stream traces just below {-y direction) the model centerline, as shown in the x-y view, 
roll up into the counter-clockwise vortex (y-z view). In contrast, the remaining lower streamlines in the x-y view 
rollup into the stronger clockwise rotating vortex. The clockwise rotating vortex is identifiable on the left side in the 
y-z view. Note that the presence of the sting is expected to influence the trajectory and strength of the counter- 
rotating vortex pair. As a result, caution must be exercised when extrapolating these observations to the wake of the 
MSL flight vehicle. Unfortunately, the presence of the sting also blocked the projection of the laser sheet through 
the lower counter-clockwise vortex and therefore prevented study of the NO entrainment into this predicted vortex. 
These observations, however, are relevant for model validation purposes since the sting was included in the 
simulations. 
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Fig. 6 3D comparison of the separated wake. PLIF data is shown as iso-surface. Stream traces predicted by 
OVERFLOW. Red arrow indicates location of dip in separated profile as described in text. 

B. RCS Jet-Flow Interaction 

The interaction of the RCS jets with the capsule wake is of interest both for aerodynamic and 
aerothermodynamic considerations. Figure 7 shows RCS B2 jet visualization with the MSL model oriented at -20° 
AoA. The image in Fig. 7a was produced using the ViDI technique. Similar to the wake-seeding technique, the RCS 
jets are seeded with NO and the laser sheet is scanned in the spanwise direction. As a result, the three-dimensional 
structure and overall trajectory of the jet is visualized. Figure 7b shows the same data visualized as an iso-surface of 
the PLIF signal (Sf= 20) using Tecplot. Note that values of Sf are in arbitrary units. The experimental PLIF image 
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shows evidence of turbulent structures on the edges of the jet’s shear layer. The scale and distribution of turbulent 
structures might be misleading, however. Because the three-dimensional jet was created by combining a series of 
uncorrelated single-shot PLIF images, large variations in spatial structure between neighboring slices can result in 
an artificially corrugated surface. Through smoothing of the experimental data in Tecplot, it is possible to filter the 
artificial small-scale structures but retain real large-scale fluid dynamic structures. Figure 6c shows the data set after 
it has been filtered with 400 passes of Tecplot’ s data smoothing algorithm. Each pass of smoothing moves the signal 
value at a node closer to the average value of the neighboring nodes. The default coefficient of 0.5 and fixed 
boundary conditions were specified in the smoothing algorithm. 



Fig. 7 Visualization of RCS Jet B2. a) ViDI visualization reproduced from Ref. [4]; b) Tecplot visualization of 
PLIF signal iso-contour surface for S f = 20 (un-smoothed); c) Tecplot visualization of PLIF signal iso-contour 

surface for S f = 10 (smoothed). 


Figure 8 compares a computed RCS jet structure with both smoothed and un-smoothed experimental PLIF iso- 
surfaces. Note that the sting, present in the experiments and computations, has been removed from the y-z view for 
clarity. The computed PLIF iso-surface is based on an empirical CFI model developed by Inman for an under- 
expanded jet [19]. The PLIF signal, 5), is calculated from the following expression: 

s f « ^,{(- 0 . 0018 )T~ 2A2 P 2 +( 46 . 0 )T~ L9 P} (1) 

In the expression, / N0 is the NO mole fraction, T is the gas temperature in Kelvin, and P is the gas pressure in 
Pascals. Since pure N 2 was used as the RCS jet fluid in the CFD simulations, a mole fraction ratio of/ NO /^N 2 = 0.05 
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was assumed. The empirical model is valid over the temperature range of 100 K to 500 K and a pressure range of 2 
Pa to 30 kPa. Rapid repopulation of the ground state, negligible saturation of the fluorescence transitions, a triple - 
Gaussian spectral overlap integral representation, and negligible quenching from 0 2 are assumed. Full details of the 
model can be found in Refs. [19] and [20]. The model was also used to perform a similar comparison between 
experiment and CFD where RCS jet interaction of the Crew Exploration Vehicle was investigated [21]. 

The computed jet structure in Fig. 8 corresponds to an iso-contour value of Sf= 0.005. For the smoothed and un- 
smoothed experimental data, the structure corresponds to a PLIF signal iso -value of Sf= 10 and S/= 20, respectively. 
The iso-value was specified for each case so that the jet structure near the nozzle exit was approximately the same 
size. Qualitatively, the shapes and trajectories of the jets predicted by OVERFLOW agree well with the experiment 
in each of the three views. Small scale turbulence structure cannot be visualized from the RANS-type numerical 
simulations as the entire turbulence spectrum was modeled. Even with data smoothing, the lack of measurements at 
some planes caused relatively large distortions in the overall shape (see y-z view). Although the jet nozzles are 
conical, both the experiment and computations show larger growth of the jet in the z- direction compared to the y- 
direction. Because the individual RCS jets are closely spaced, it is likely that the jet-jet interactions promote greater 
expansion in the z- direction compared to the y- direction. The initial jet trajectory direction is away from the sting (- 
y) but then reverses direction back towards the sting further downstream. The redirection of the jet towards the sting 
occurs from the interaction of the jet with the high-speed outer flow. 
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Fig. 8: Comparison of computed S f = 0.005 (first row), smoothed experimental S f = 10 (second row) and un- 
smoothed experimental S f = 20 (third row) RCS Jet iso-contour surfaces. 


Typically, visualization of the jet structure from CFD involves plotting an iso-surface of a single flow parameter 

such as pressure, Mach number, or mixing fraction. For example, Fig. 9a shows an iso-surface of/ NO = 0.038 (mole 

fraction of NO is a parameter equivalent to jet fluid mixing fraction for this set of computations). By combining 
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multiple parameters using Eq. (1), CFI creates an alternative visualization of the jet structure by generating iso- 
surfaces of computed fluorescence intensity. Figure 9b shows the influence of the CFI model on the visualized shape 
of the jet structure with a selected iso-surface of Sf= 0.005. The/ N0 iso-surface value was selected to match the jet- 
structure and size of the CFI model iso-surface near the nozzle exit. Note that although the shape and trajectory of 
the jet structures are similar, the CFI model predicts an iso-surface that persists further downstream of the capsule. 



Fig. 9 Effect on RCS jet structure of plotting a CFD flow parameter versus using a CFI model: a) Iso -surface 

of Xno = 0.038; b) Iso-surface of S f = 0.005. 


Upon further analysis of the simulation solution, it was found that shock compression of the RCS jet in the 
capsule wake results in an increase in gas density and ultimately an increase in PFIF signal intensity. The effects of 
gas density are indirectly accounted for in Eq. (1) from the combination of pressure and temperature terms. Figure 
10 shows the interaction of the RCS jet with the outer-flow, interaction with the first oblique shock wave that 
originates from the sting, and interaction with the second oblique shock that originates from the wedge further 
downstream on the sting. As shown in the figure, each of these interactions results in an increase in PFIF signal for a 
constant NO mole fraction. Prior to the shock wave interactions the jet structure appears to be mostly dominated by 
the NO mole fraction distribution as it forms in the subsonic wake near the capsule surface (see Fig. 9). 
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Fig. 10 CFI result showing dependence of fluorescence signal on factors other than mole fraction of NO. Iso- 
contour of S f = 0.005 colored by contours of Sfy no . 


Figure 11 shows the jet trajectory for both the simulation and experiment based on post-processing of the 
respective PLIF signals. Using a Tecplot macro script, each point is calculated from signal-weighted average over a 
y-z plane at a given x position: 


IX,* IX,,.-, 
' IX, • * ' IX, 


In Eq. (2), y c and z c are the y and z coordinates corresponding to the center of the jet. The trajectory is constructed as 
the calculation is performed at subsequent downstream y-z planes. A similar technique was used to calculate the 
two-dimensional trajectory from two dimensional PLIF images of the RCS jets for an Apollo -geometry capsule [22]. 
Figure 12 shows a direct comparison between the x-y and x-z jet trajectories based on the experimental and 
computed flow visualization data. In the experimental flow visualization data, a threshold of Sf = 20 was applied to 
Eq. (1) to remove artificial effects of noise on the calculation far away from the jet center. 


18 


American Institute of Aeronautics and Astronautics 




Fig. 11 RCS jet trajectory based on computed (left) and experimental (right) signal iso-contour surfaces. 


Also shown in Fig. 12 is the effect of applying a threshold to the calculation of the trajectory. Results indicate 
that in the simulation, the trajectory near the nozzle exit is unaffected by the threshold level. There is, however, 
some sensitivity of the calculation to threshold levels further downstream. As the NO diffuses away from the jet core 
and peak concentration levels decrease, lower thresholds are required to accurately predict the jet trajectory. Overall, 
the agreement between the simulation and experiment is marginal. Errors are most likely caused by the threshold 
method applied to the experimental data and some inaccuracy of the empirical PLIF model. Lower thresholds could 
not be applied to the experimental data since the jet is not centered on the measurement volume and increased noise 
levels near the sting bias the trajectory in the +y direction. Better signal-to-noise measurements in the experiment 
would improve the accuracy of the trajectory calculation. 
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Fig. 12 Comparison of the computed and experimental jet trajectory shown in Fig. 11. The effect of threshold 

levels on the predicted trajectory is shown. 


C. Velocity Measurement Comparison 

NO PLIF molecular tagging velocimetry was used to characterize the streamwise velocity field in the wake flow 
with and without the RCS B2 jet pair activated. When the RCS jet pair was on, NO was seeded both into the RCS 
jets and from the junction between the model and sting. When the RCS jet was off, NO was seeded only from the 
junction between the model and sting. For each RCS jet operational mode (jet on or jet off), the laser sheet traversed 
a set of six equally spaced spanwise distances, stopping at each location for approximately 2 seconds to allow for 20 
images to be captured at each location for statistical analysis. Through post-processing, the quantitative distribution 
of streamwise velocity along each of the MTV lines was measured. A contour map of streamwise velocity was 
created through the interpolation of each of the MTV lines onto a two dimensional plane using Tecplot. A sample of 
one of the contour maps of streamwise velocity (V x ) without the RCS jets activated is shown in Fig. 13 (right) along 
with the corresponding prediction from OVERFLOW (left). In this figure, the laser sheet is positioned at y = -39 
mm, which is the spanwise position measured from the model centerline. This spanwise position is roughly one 
nozzle diameter away from the activated RCS jet pair. 
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Since the RCS jet is not activated, the magnitude of streamwise velocity in the separated wake region is very low 
compared to the speed of the supersonic outerflow. Furthermore, NO gas seeded through the sting is unable to 
diffuse through the shear layer into the high-speed outerflow. As a result, the field of view of the measurement is 
limited to the size of the wake. With low signal-to-noise levels associated with the wake-seeding technique and a 
limited number of data points obtained during each run, the measurement error is of the same order of magnitude as 
the measurement itself. Although the majority of velocity measurements were low in magnitude, higher velocity 
magnitudes were measured in the top shear layer. Unsteadiness in the wake region, which is confirmed from the 
visualization images in Fig. 6, will cause an error in the position of the shear layer. MTV measurements will over- 
predict the position of the shear layer since only high signal (low magnitude velocity) measurements will be 
included in the average measurement. NO velocity measurements in the shear layer and in the outerflow could be 
obtained if the NO gas were seeded upstream from the capsule’s forebody. The problem, unfortunately, with 
forebody seeding is that there is a risk that the seeded gas will inadvertently perturb the attached boundary layer 
flow. Since the velocity distribution in the wake is relatively uniform, only the overall velocity magnitude and shear 
layer location can be compared to the experiment. Further comparisons between the jet-off CFD and experimental 
data in four additional spanwise planes are available in the literature [23]. 



Fig. 13 Contour map of CFD predicted (left) and experimentally measured (right) streamwise velocity ( V x ). 
Laser sheet positioned a ty = -39 mm from the model centerline. RCS jet is deactivated. 
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Figure 14 shows a contour map of the computed (left) and experimentally measured (right) streamwise velocity 
at the same spanwise position as that shown in Fig. 13 but with the RCS jets activated. Unlike the jet-off case, large 
magnitudes of streamwise velocity are measured in the wake. Because the RCS jets are internally seeded with NO 
gas, good signal -to-noise levels are achieved within the jet. As a result, a comparison between the computed and 
experimental RCS jet structure can be made. The simulation under-predicts the streamwise velocity within the jet 
core by about 10% but has good agreement on the overall shape and size of the jet. High magnitudes of velocity on 
the model surface are not expected and are an artifact of the laser scatter in the experiment. 



Fig. 14 Contour map of CFD predicted (left) and experimentally measured (right) streamwise velocity ( V x ). 
Laser sheet positioned aty = -39 mm from the model centerline. RCS jet is activated. 


The simulation does predict a region of high-velocity gas above the main jet structure, which although subtle, 
also appears in the experimental data. Further analysis of the CFD simulation reveals that this high-velocity gas 
region corresponds to the core of the counter-clockwise rotating vortex described in the Wake Flow Visualization 
section. It is possible that the interaction between the low-pressure vortex core and the RCS jet leads to unwanted 
changes in the surface pressure distribution with regards to control authority. This interaction is very sensitive to the 
model’s AoA as both the strength and location of the counter rotating vortex are affected. Further two-dimensional 
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comparisons between the jet-on CFD and experimental data in three additional spanwise planes are also available in 
the literature [23]. 

Figure 15 shows the same computed velocity field as the left figure in Fig. 14 but with dashed lines indicating 
the locations of the individual measured velocity profiles. The individually measured streamwise velocity profiles 
along these lines corresponding to the jet-off and jet-on cases are shown in Fig. 16 and 17, respectively, along with 
the CFD predictions. The width of the symbols in Figs. 16 and 17 indicates the 95% confidence uncertainty bounds 
of the measured mean velocity at that location. For the most part, the predictions are within the uncertainty bands of 
the experimentally measured velocities. 



Fig. 15 Contour map of computed streamwise velocity ( V x ) indicating the positions of the measured velocity 

profiles. 


Confirming what was observed qualitatively in Fig. 13, comparisons between computed and experimental 


velocity profiles show that the simulations under-predict the height of the shear layer location. Placement of the 


experimental MTV data relative to the model was based on the location in the images of laser scatter off the surface 


of the model. As described above it is likely that there is a biased error in the location of the shear layer because of 
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the unsteadiness in the wake structure and the systematic rejection of low-signal measurements when the outerflow 
penetrates the region of interest. 



Fig. 16 Vertical (z) distribution of streamwise velocity (K x ) for the conditions of Fig. 13 at y = -39mm. RCS jet 
is deactivated. Width of bars indicates uncertainty in the mean values. (CFD-solid line; MTV-points) 


As expected the measurement error within the RCS jet core is relatively low because of the high signal -to-noise 
ratio in the PLIF images. There is good agreement between the predicted and experimental vertical location 
corresponding to peak velocity. Confidence in the vertical placement of the MTV data indicates that the poor 
prediction in the shear layer location is likely due to a bias error in the MTV measurement. It is also possible that the 
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prediction accuracy of the separation location is limited by the RANS-type turbulence modeling. RANS-type 
turbulence models traditionally have difficulty predicting flows with streamline curvature and flow separation [9]. 
Further analysis is needed to determine the discrepancy in the shear layer location. Another discrepancy between the 
simulation and experiment is the prediction of the peak magnitude of velocity. At each streamwise location, the 
simulations under-predict the magnitude of peak velocity. Although additional comparisons in the literature indicate 
that the peak velocity is sensitive to the spanwise location [23], the peak velocity extracted from the MTV 
measurements is always larger than the predicted peak values. The under-prediction of peak streamwise velocity 
magnitudes can be partially attributed to an error in the predicted jet orientation. For example, the initial x-z and x-y 
trajectory measurements (75 mm < x < 95 mm) shown in Fig. 12 indicate that the predicted ratios of streamwise-to- 
vertical (V z /V x ) and streamwise-to-spanwise (Vy/V x ) velocity components are lower than the ratios calculated from 
the experiment. As a result, the predicted total velocity magnitude could be accurate but include errors in the 
predicted velocity component ratios (i.e. VJV X and Vy/V x ). The observation in the trajectory described above occurs 
over a streamwise range (75 mm < x < 95 mm) that is coincident with measurements shown in Fig. 17. It is also 
possible that an NO LIF quenching gradient exists in the streamwise direction and caused a bias error in the 
calculation of the velocity magnitudes from the experimental MTV data [16]. Furthermore, the MTV velocity 
calculations might be sensitive to the maximum signal intensity threshold criteria used to reject large signals in the 
first camera gate that artificially deplete signals in the second camera gate. Depletion of signals in the second camera 
gate will bias the velocity measurement to higher values. 
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x = 78.7 mm 



Streamwise Velocity (V x , m/s) 


Fig. 17 Vertical ( z ) distribution of streamwise velocity (F x ) for the conditions of Fig. 14 at y = -39mm. RCS jet 
is activated. Width of bars indicates uncertainty in the mean values. (CFD-solid line; MTV-points) 


VI. Discussion 

There are some limitations in the experimental PLIF technique and CFD modeling that have introduced errors 

into the comparisons of the MSL wake and RCS jets. Experimentally, the unsteadiness of the capsule wake 

introduces bias errors into the position of the shear layer and in the size of the recirculation zone. Although there is a 

risk of perturbing the flow, NO seeding from the capsule forebody could improve the visualization and MTV 

measurements of these flow structures. The MSL model described in Ref. [4] has ports on the forebody to enable 
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this type of seeding. The MTV velocity measurements could also be improved if the NO wake seeding flow rate 
were increased above levels reported in Ref. [4]. A higher concentration of NO in the wake would improve the 
signal-to-noise ratio in the tagged lines. In tests where the RCS jets are activated, the MTV technique could be 
improved if the NO concentration in the RCS jets were decreased to avoid saturation of the camera’s CCD. Several 
improvements could be made to the numerical modeling. With an increase in computational resources, it would be 
possible to run the simulations in an unsteady mode and possibly use a more sophisticated turbulence modeling 
approach such as LES. With unsteady numerical simulation data, a comparison to the single -shot experimental PLIF 
data could be made. Finally, the CFD predictions could be improved if a more sophisticated CFI model were used. 

The simulation solutions have identified important aerodynamic structures in the wake of the MSF vehicle that 
are not currently captured by the experimental PFIF dataset. For example, a counter-rotating vortex pair that entrains 
NO from the sting and transports it far downstream was identified. Since the laser sheet was orientated parallel to 
the flow and was positioned in the near wake, these structures were largely undetected experimentally. Assuming a 
negligible interaction with the model sting, the circulation of these trailing vortices could be related to the 
aerodynamic forces acting on the vehicle through a process analogous to Prandtl’s lifting line theory. The circulation 
of these trailing vortices could be extracted experimentally in a carefully designed PFIF experiment. By re- 
orientating the laser sheet and cameras in the streamwise direction, the x-component of vorticity can be measured. A 
two-component (cross -pattern) MTV velocity measurement would capture vorticity in the separation regions along 
the centerline. This method could provide an alternative, non-intrusive diagnostic to measure model forces. 

VII. Conclusions 

A validation study of the OVERFFOW code has been performed through qualitative and quantitative 
comparisons of the predicted Mars Science Faboratory wake and reaction control system jet flow field to an 
experimental planar laser-induced fluorescence dataset. Overall there is good agreement between the experiment and 
simulation with regards to the size and shape of the separated recirculation zones in the capsule wake. A counter- 
rotating vortex pair was identified from the numerical simulations, the effect of which was not deduced from the 
experiments alone. The simulations predict entrainment of nitric oxide seeded from the sting and jets into the vortex 
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cores, allowing for potential experimental fluorescence measurements to be obtained in the far field. A 
computational flow imaging model was applied to the numerical solution to make comparisons between the 
predicted and experimental jet structure. The effect of smoothing on the experimental data was investigated and the 
processed images were compared to the simulation from several different views. The simulation solution was used 
to investigate the limits of the fluorescence signal model. It was found that near the capsule surface the jet structure 
based on the fluorescence signal is dominated by mole fraction of nitric oxide. Further downstream the interaction of 
the jet with a series of shock waves influences the signal levels and increases the perceived jet size and length 
compared to structures based on iso-surfaces of nitric oxide mole fraction. For the first time, a three-dimensional jet 
trajectory was extracted from both the experimental and computational flow images and compared. Agreement 
between the experiment and simulation based on the jet trajectory was marginal. Experimental molecular tagging 
velocimetry measurements were compared both qualitatively and qualitatively for jet-on and jet-off cases at several 
spanwise positions in the capsule wake. There was good agreement in the regions with high signal-to-noise ratio, 
such as in the jets. Low signal-to-noise levels in the experimental fluorescence data away from the jets restricted the 
comparison for the jet-off case. Distributions of high-velocity gas predicted by the simulation and confirmed in the 
molecular tagging velocimetry measurements were found to coincide with the core of the streamwise trailing vortex. 
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